1 



EXTRACTING STELLAR PARAMETERS FROM OBSERVED SPECTRA: THE ROLE OF 
CROSS-CORRELATION AND MINIMUM DISTANCE METHODS 
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SCI-4, Thevenin et al. 2003, Katz et al. 2004). Note that 
these approaches are relevant at the end of the GAIA mis- 
sion, but at the start of the mission virtually nothing will 
be known about the values of stellar parameters, so that 
one would not know which kind of template to use for a 
radial velocity determination. Also, the wavelength cov- 
erage of GAIA spectra is much smaller than that of a typi- 
cal Echelle (ELODIE) type spectrum. So we attempt here 
to determine both the best radial velocity and the values 
of parameters of the stellar atmosphere simultaneously. 
Data from the RAVE experiment, which are very similar 
in both wavelength coverage and resolution to those of 
GAIA's RVS, will be used to illustrate the point. Note 
that since virtually nothing is known about the RAVE tar- 
gets it is obvious to use such an approach. 



One may consider cross-correlation routines applied to 
stellar spectra as a way to determine not only the radial 
velocities but also to assess the properties of stellar atmo- 
spheres. We show that this is not the case. High value 
of the cross-correlation coefficient between the observed 
and a template spectrum does not mean that the template 
corresponds to correct values of stellar temperature, grav- 
ity or metallicity. Cross-correlation operates in normal- 
ized space, so it is not sensitive to scaling or shifting of 
stellar flux. It is argued that cross-correlation is useful 
to efficiently determine the velocity shift between the ob- 
served and the template spectrum. But minimal distance 
methods need to be used to determine which template 
gives the best match to the observed spectrum, and so de- 
rive values of physical parameters. To illustrate the point 
we use data from the RAVE project which are being ob- 
tained at very similar spectral resolution and wavelength 
domain as will be the case for GAIA's RVS. 

Key words: Stars: fundamental parameters — Methods: 
data analysis. 



1. INTRODUCTION 



In the attempt to extract values of physical parameters 
from stellar spectra one first tries to determine the value 
of radial velocity. If the spectra have a wide wavelength 
coverage, like in the case of ELODIE spectra, it can be 
assumed that a fixed spectral mask gives radial veloci- 
ties with sufficient accuracy (Katz et al. 1998). Analy- 
sis of GAIA radial velocity spectrometer (RVS) data will 
profit from quite accurate knowledge of stellar param- 
eters obtained from on-board photometric observations. 
So all simulations of GAIA radial velocity determination 
assumed that the choice of a proper template is given by 
complimentary photometric data (Katz et al. 2004, Mu- 
nari et al. 2003, Zwitter 2002). Alternatively it was as- 
sumed, while assessing the capabilities of GAIA spectra 
to independently determine the values of stellar parame- 
ters, that the radial velocity is a-priori known (ES A-2000- 



2. OBSERVED SPECTRA 



RAVE (RAdial Velocity Experiment) project (Steinmetz 
2003, Munari et al. 2004) is an ambitious interna- 
tional collaboration aimed at conducting an all-sky spec- 
troscopic survey of Galactic stars. The participating 
countries are Australia, Canada, France, Germany, Italy, 
Japan, Netherlands, Slovenia, Switzerland, UK, and 
USA. The spectra allow us to measure radial velocities 
as well as metallicites and other stellar parameters for a 
large number of stars using the 1 .2-m UK Schmidt tele- 
scope at the Anglo-Australian Observatory. Data col- 
lection for this project started in April 2003 with more 
than 40,000 spectra observed to date. Spectral resolv- 
ing power and wavelength range (8413-8746 A) are very 
similar to those of the GAIA Radial Velocity Spectro- 
graph. 



3. LIBRARY OF SYNTHETIC SPECTRA 



Observed or synthetic spectra can be used for classifica- 
tion. The advantage of the former is that they are "real", 
so avoiding any simplifying assumptions used in syn- 
thetic spectra calculation (Munari et al. 2001). The lat- 
ter have the advantage of a more uniform coverage of 
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Figure 1. Fits of templates with the best value of correlation coefficient R are not necessarily those with the smallest value 
of the standard deviation { 0-C) between the observed and template spectra. The later is given in units of continuum. Each 
point presents results of a cross-correlation and standard deviation for one template from the database. Symbol type marks 
the difference between the calculated value of radial velocity for a given template and the one for a minimal 0-C template: 
less than 0.2 km (•}, up to 2 km (o), or more (A). 



the parameter space. At present we are using a subset 
of a synthetic Hbrary of nearly 60.000 spectra with suit- 
able resolving power (Zwitter, CastelU & Munari 2004) 
which were calculated from the latest generation of Ku- 
rucz models. The grid is characterized by the following 
ranges of stellar parameters: 3500 K < Toff < 47500 K, 
0.0 < log5 < 5.0, -3.0 < [M/H] < +0.5, < Vrot < 
500kms-\ ^ = 2kms-^ 



4. USING CORRELATION ROUTINES ONLY 



Correlation routines, like the IRAF's rvsao package 
(Kurtz & Mink 1998) are a well estabUshed way to cal- 
culate radial velocity using a fixed stellar or galaxy spec- 
trum template. Their figure of merit is based on the value 
of correlation coefficient R (Tonry & Davis 1979). But 
one should not conclude that the template spectrum with 
the highest value of R also gives the best fit to the ob- 
served spectrum. Fig. 1 shows results of analysis of a 
single observed spectrum using a set of template spectra. 
The maximal value of the correlation coefficient R oc- 



curs at the "stem" of area covering a caravel ship shape, 
though the minimum of O — C occurs at its bottom. The 
reason is that the correlation is defined using normalized 
quantities. Therefore it is not sensitive to any kind of con- 
stant shift or spread (of the kind x -\- a and x*a where a 
is a constant value or bias). This explains why it is only 
sensitive to high variations in the line profile and is not 
suitable for minimization. 

We can also say that in cross-correlation one tries to find 
the best alignment of the lines, in a way irrespective of 
their depth. This is strictly valid for isolated single lines. 
The effective wavelength of a blend of nearby lines can 
shift sideways and so worsens the cross-correlation only 
with a large change of the temperature and gravity. 

Figure 2 plots the observed spectrum which was used in 
Figure 1 and one of the templates from the library. The 
match is poor even though the template shown features 
the highest value of the correlation coefficient R among 
all templates in the library (the point is at the tip of the 
stem of the caravel shape in Figure 1). Cross-correlation 
alone is clearly not enough to find the template with the 
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Figure 2. Match between the observed (thin-line) and template spectrum (thick-line) is poor, even though this is the 
template which yielded the highest value of the correlation coefficient R in Figure 1. 



best match. 



6. RELEVANCE OF RAVE, GAIA AND OTHER 
AUTOMATED CLASSIFICATION PROJECTS 



5. ADDING A MINIMAL DISTANCE METHOD 



The situation can be improved by supplementing correla- 
tion analysis with a proper minimization scheme (Figure 
3). The reason is that normalized squared distance is sen- 
sitive to line depth and shape as well as continuum level, 
and so to the values of stellar parameters. 

In our case we first use a correlation routine to calculate 
the velocity shift between the observed spectrum and a 
template from the synthetic library. But the goodness of 
fit is not based on the value of the correlation coefficient 
R. Instead we calculate the standard deviation between 
the observed spectrum and the velocity- shifted template. 
The template which yields the smallest standard deviation 
is assumed to give a fair representation of parameters of 
the stellar atmosphere. 



Values of several stellar parameters are primarily based 
on spectral analysis. These include metallicity, abun- 
dances of individual elements and stellar rotation veloc- 
ity. The values of effective temperature and surface grav- 
ity can also be derived, or in the case of GAIA checked, 
by spectral analysis. Both GAIA and RAVE projects 
will provide too many spectra to be analyzed by hand. 
Automated procedures, such as the one proposed here, 
yield useful results. Any other information available, in 
the case of GAIA these would be results of photometric 
and astrometric observations, can be seen as an additional 
constraint on the relevant range of stellar templates to be 
considered for a given star. If an acceptable match could 
not be found despite a high S/N of the observed spectrum, 
this imphes that any kind of pecuUarity is present in the 
observed object, or that the observed star is a double lined 
spectroscopic binary. 
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Figure 3. The match between the observed (thin-line) and template spectrum (thick-line) is much better when a suitable 
minimization scheme is used as a figure of merit. The template and observed spectrum correspond to the lowest O — C 
point in Figure 1. 



7. CONCLUSION 



Cross correlation is sensitive to position of the lines and 
so radial velocity. It is however not sensitive to line depth, 
so it is not to be used to determine the values of stellar at- 
mosphere parameters unless blends of nearby lines are 
driving the correlation. The situation gets worse with low 
S/N spectra and when sharp hues (e.g. Ca II IR triplet) are 
present. By adding a suitable minimum distance scheme 
which is sensitive to equivalent widths (Thevenin & Foy 
1983) the values of both radial velocity and stellar pa- 
rameters like temperature, gravity, metallicity and rota- 
tion velocity can be determined reliably. 
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